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Abstract. We study Kerr nonlinear resonators (KNR) driven by a continuous wave field in quantum 
regimes where strong Kerr interactions give rise to selective resonant excitations of oscillatory modes. 
We use an exact quantum theory of KNR in the framework of the Fokker-Planck equation without any 
quantum state truncation or perturbation procedure. This approach allows non-perturbative consideration 
of KNR for various quantum operational regimes including cascaded processes between oscillatory states. 
We focus on understanding of multi-photon non-resonant and selective resonant excitations of introcavity 
mode depending on the detuning, the amplitude of the driving field and the strength of nonlinearity. The 
analysis is provided on the base of photon number distributions, the photon-number correlation function 
and the Wigner function. 

PACS. 42.50.Ar - 42.50.Hz 


1 Introduction 


It is known that a cavity mode strongly coupled to a non¬ 
linear element leading to direct photon-photon interac¬ 
tions behaves like a one-mode oscillator with an anhar- 
monic energy ladder. Large nonlinearities make the pho¬ 
ton energy spectrum non-equidistant and well resolved, 
therefore spectroscopic selective excitations of transitions 
between Fock states become possible. Currently existing 
and recently proposed oscillatory systems realizing photon- 
photon interactions contain a large number of various non- 
linearities including Kerr-nonlinearity. A Kerr nonlinear 
resonator (KNR) can serve as an example. Selective exci¬ 
tations of KNR can be realized by tuning the frequency of 
the driving field to the sum of the resonator frequency 
and the Kerr nonlinear shift of the energy levels. The 
anharmonic, Jaynes-Cummings ladder has been recently 
probed in atomic and superconducting cavity quan¬ 
tum electrodynamics (cQED) systems. It has been ob¬ 
served that non-classical correlations between photons be¬ 
ing transmitted through the cavity can result from such 
anharmonicities, leading to the fundamental phenomena 
of photon blockade (PB) and photon-induced tunneling. In 
PB, the capture of a single photon by the system affects 
the probability of a second photon being admitted. PB 
has been proposed for many applications, including the 
generation of single photons, quantum communications, 
where strong antibunching of photons and sub-Poissonian 
statistics is required, high precision sensing and quantum 
metrology. 


Originally, PB was proposed in a Kerr-type nonlin¬ 
ear cavity via giant photon-photon interactions and 
was first observed in an optical cavity coupled to a single 
trapped atom [^. PB has also been predicted in cQED [^, 
and recently in circuit QED with a single superconduct¬ 
ing artificial atom coupled to a microwave transmission 
line resonator IgII^. Going further, PB was experimentally 
demonstrated for a photonic crystal cavity with a strongly 
coupled quantum dot [s], and was predicted in quantum 
optomechanical system ^9|IQ| . An analogous phenomenon 
of phonon blockade was predicted for an artificial super¬ 
conducting atom coupled to a nanomechanical resonator 
II , and the polariton blockade effect due to polariton- 
polariton interactions has also been considered 12 . Re¬ 


cently, one-photon PB was considered for dispersive qubit- 
field interactions in a superconductive coplanar waveguide 


cavity 


with time-modulated input [14], and for the 


train of pulses . Two-photon blockades as well have 


been investigated 15 17 . It is commonly believed that 


photon blockade necessarily requires a strong Kerr non¬ 
linearity for a photonic mode, magnitude of which should 
well exceed the mode decay rate. It has been also shown 
that photon blockade could be achieved in nanostructured 
cavities with third-order nonlinear susceptibility 18 . An 


alternative unconventional photon-blockade effect was re¬ 
cently proposed with a weaker Kerr nonlinearity than its 
conventional counterpart in a system consisting of two 
coupled photonic cavities or from suitably engineered cou¬ 
pled modes in one cavity 19 . The strong photon-photon 


correlation in this scheme is attributed by the destructive 
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quantum interference between distinct driven-dissipative 
pathways 21 . Based on this mechanism, many dif¬ 


ferent systems are proposed to realize the unconventional 
photon-blockade [2^29 , moreover three oscillatory mode 


systems have also been considered 30 


used to reach strong nonlinear regimes. Particularly, the 
regimes of weak x /7 << 1 , strong x /7 > 1 and very 
strong x /7 >> 1 nonlinearities have been demonstrated 
for a Josephson junction embedded in a transmission-line 


resonator by adjusting the parameters of the circuit 44 


In all cited papers, the regimes of PB in cavity and 
determination of optimal values of the detuning and the 
nonlinear coupling constants are analyzed in a simplified 
approach based on the wave function study. Here, mode 
losses are treated at the non-hermitian hamiltonian level 
and the wave function expands on a Fock-state basis. As¬ 
suming weak optical pumping the dynamics of the lowest 
Fock state is only analyzed in a steady state. In this way, 
the master-equation solutions are obtained through ma¬ 
trix elements of Fock states analytically as well as numer¬ 
ically for the lowest states in steady state regimes. This 
approximation of Fock state truncation in a steady-state 
allows to consider a Kerr-type nonlinear resonator or a 
system of coupled resonators for the case of sufficiently 
weak driving external fields. In this regime the photon- 
number distributions are also discussed [^,32 . Besides, 
the processes of re-excitation and cascaded decay of os¬ 
cillatory states are neglected in this approximation. This 
approach as well as the nonlinear quantum scissor effects 


46 


were discussed earlier |^(see, also the review paper 34 


Photon-number effects, multi- photon blockade and gen¬ 
eration of Fock states in KNR have been studied numeri¬ 
cally in the full time-evolution of the master eq uation o n 
the base of the quantum state diffusion method 


The paper is arranged as follows. In Sect. we present 
the quantities of interest for KNR by using the poten¬ 
tial solution of the Fokker-Planck equation in complex P- 
representation. In Sect. we analyze multi-photon excita¬ 
tions in KNR on the base of the mean photon number, the 
photon-number state populations, the second-order pho¬ 
ton correlation function and the Wigner functions. This 
analysis includes a study of non-resonant cascaded exci¬ 
tations as well as selective resonant excitations of KNR 
depending on the detuning, the amplitude of driving field 
and the strength of the Kerr nonlinearity. We summarize 
our results in Sect. [U 


2 Exact quantum treatment of KNR: short 
description 

The Hamiltonian of Kerr-nonlinear resonator driven by a 
monochromatic field in the rotating-wave approximation 
reads as: 


35 38 


In this paper, we analyze selective excitations of photon- 
number states in KNR driven by a continuous wave field 
in a steady state regime for strong Kerr nonlinearities. 

The novelty here is that we analyze selective excitations 
in KNR using exact quantum treatment of driven dissi¬ 
pative KNR in the framework of the Fokker-Planck equa¬ 
tion in complex P-representation 39 -1^ . This approach 
gives a full quantum description oTthe Kerr dissipative 
resonator for quantum noises of arbitrary strength with¬ 
out any state truncation procedure and hence allows us 
to find fundamental limits of PB due to dissipative and 
decoherence effects. Moreover, the approach allows non- 
perturbative consideration of KNR for various quantum 
operational ranges including cascaded processes between 
oscillatory states and multi-photon selective resonant ex¬ 
citations. The last case is realized by tuning the frequency 
of the driving field for wide ranges of the parameters. We 
clarify the selective excitations of the oscillatory mode 
leading to multi-photon blockades considering photon-numb efollowing form: 
effects on the base of the mean photon number, the prob¬ 
ability distributions of photon numbers, the second-order 
correlation function of photons and the Wigner function 
of the cavity mode. 


Hq = TiAo) a + hxci^^ci^ + h{Oa^ + i?*a). 


( 1 ) 


where a'*', a are the oscillatory creation and annihila¬ 
tion operators, x is the nonlinearity strength and A = 
cjo ~ ^ is the detuning between the oscillatory frequency 
and the frequency of the driving field. The Hamiltonian 
Q describing an anharmonic oscillator has been proposed 
for wide range of physical systems, including optical fibers, 
nano-mechanical oscillators, various Josephson junction 
based devices, quantum dots, quantum scissors, etc. 

The full Hamiltonian of the system reads as = 
Hq + ^loss’ term J^Iqss ~ + a+T is re¬ 

sponsible for linear losses of the oscillatory mode, due to 
couplings with the heat reservoir operators giving rise to 
the damping rate 7 . We trace out the reservoir degrees 
of freedom in the Born-Markov limit assuming that the 
system and the environment are uncorrelated at initial 
time t = 0. This procedure leads to a master equation for 
the reduced density matrix of the oscillatory mode in the 


dt 


^[Ho,p]+ ^ LipL+-htLip-^pLtLi. ( 2 ) 

i=l,2 


The important parameter responsible for KNR is the 
ratio x /7 between the parameter of the Kerr-type non¬ 
linearity and the damping of the oscillatory mode (the 
photon decay rate). As a rule, the efficiency of quantum 
nonlinear effects requires a high nonlinearity with respect 
to the dissipation, posing a severe technological challenge. 
Nevertheless, resonators with Josephson junctions can be 


Here Li = + l) 7 a and L 2 = ^/Nxa'^ are the 

Lindblad operators, 7 is the dissipation rate and N de¬ 
notes the mean number of quanta of the heath bath.To 
study the pure quantum effects we focus on the cases of 
very low reservoir temperatures with N approximately 
equal to zero. The master equation is then transformed 
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into a Fokker-Planck equation. In this simplest case ana¬ 
lytical results for the Kerr-nonlinear dissipative resonator 
in a steady state have been obtained in terms of the solu¬ 
tion of the Fokker-Planck equation for the quasi-probability 
distribution function P{a, a*) in the complex P-representation 
39 . This approach based on the method of potential equa¬ 
tions leads to the analytic solution for the quasi-probability 
distribution function P{a, a*) within the framework of the 
exact nonlinear treatment of quantum fluctuations. On 
the whole, the various moments of mode operators as well 


as the photon distribution function were obtained 40 42 


For instance, the photon number probability distribution 
function P^ = {n\p\n) can be expressed in terms of the 
complex P-representation as follows 


Pn = - 
n 


- f f dadQ;*e-““*Q;"a*”P(a,a*), (3) 

nj Jc 


where C is an appropriate integration contour for each 
of the variables a and a*, in individual complex planes. 
The details of calculations involving contour integrations 
have been presented in 4^42 . The results for the mean 
photon number and the photon number probability dis¬ 
tribution function read as 


_ |£rr(A)r(A*) g. 


\2k 


(n) = {a^a) = 


12 ^ 


oF2(A + 1; a* + 1; z) 


( 4 ) 


(Z\ + X)" +(7/2)2 oF2(A;A*;2|5|2|) 


( 5 ) 

where 5 = f2/x, A = (7 + iZ\)/ix, and 0 F 2 is the hy¬ 
pergeometric function 




n!oF 2 (A; A*; 2\e\^) ^ k\r{k + n + \)r{k + n + A*) ’ 


Fig. 1: Photon numbers in dependence on the detuning 
and the driving field amplitude. The parameters are as 
follow: x /7 = 2 (a) and x/t = 20 (b). 


Here Jy is the Bessel function, A = (7 + iA)lx^ ^ = 
i7/x, N is the normalized constant and the amplitude 
(a = X + is a complex c-number variable correspond¬ 
ing to the operator a. For further calculations we use the 
representation of Bessel function in terms of Gamma func¬ 
tions and the hypergeometric function 


oF2(a;6; z) = ^ 


k=0 


z'^r{a)r{b) 

k\r{k + a)r{k + b)' 


(6) 


Mz) = 


In this approach, the second-order correlation function 
of photon-numbers for zero-delay time where. 


y 2 )(o) = 


< a^a^aa > 


< a^a >2 
is calculated in the following form 


( 7 ) 


(-zY 1 


oFi{;b;z) = Y, 


fe =0 


z'^r{b) 

k\r{k + b)' 


( 10 ) 


( 11 ) 


( 2 ). . ^ |£|^r(A)r(A*)oF2(2 + A;2 + A*;2|£p) 1 

^ ^ ^ r(2 + A)r(2 + A*)oF2(A;A*;2|£|2) {ny 

(8) 

Of particular interest is the analysis of phase-space 
properties of cavity mode in the framework of the Wigner 
function derived from the Fokker-Planck equation of KNR 
in a steady-state regime. The Wigner function involving 
quantum noises with terms of all orders of the perturba¬ 
tion theory reads as 


W(a) = Vexp-2|«l' 


J\-iW-8as) 

(a*)(^-i )/2 


( 9 ) 


We then calculate the Wigner function in the following 
form 


W{a) = Arexp“2“" 


(- 4 e)"^^E 


{ 2 ae)’‘ 

r{X + k)k\ 


2 


( 12 ) 


This approach is applied for quantum operational regimes 
of KNR consisting of non-resonant and resonant excita¬ 
tions of the cavity mode by an external coherent field be¬ 
low. 
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(b) 


Fig. 2 : (Color online) Populations of photon numbers de¬ 
pending on the detuning (a); cross-sections of the plot for 
particular detuning values (b). 


3 Multi-photon excitations in a Kerr 
dissipative resonator 

In the absence of any driving, states of the nonlinear os¬ 
cillator are photon-number states \n) which are spaced in 
energy En = Eq ^ Tiuj^n + hxn{n — 1 ) with n = 0 , 1 ,.... 
The levels form an anharmonic ladder with a nonlinear 
parameter that is given by F^ 2 i — ^lo = The energy 
spectrum can be probed with response of the system to the 
driving field when the driving frequency is tuned. In this 
case, spectral lines of the system are resolved if the non¬ 
linear shifts of oscillatory energy levels are larger than the 
state line-widths, i.e. x /7 > 1 . According to the formula 
of the energy spectrum resonant frequencies nhuOn = Eno 
of the multiphoton transitions | 0 ) ^ |n) are derived in 
the following form ujn = ojq x(n — 1). Thus, in order to 
excite a state with one photon {Eio = Tiujq) the driving 
frequency is fixed to the resonant frequency uJi = For 
the two-photon transition | 0 ) ^ | 2 ), E 20 = 2?iujq + 2 %, 
and the resonant frequency is Cc ;2 = ^0 + A, while for the 
three-photon transition |0) ^ |3), = SficJo + 6 %, the 

resonant frequency is = ujq E 2 %. Hence, the detun¬ 
ing values at the resonant frequencies can be written as 
An = uJo — uJn = ~ !)• This analysis is correct if the 

driving field amplitude is weak. Otherwise the spectrum 
of mode excitations will contains frequencies correspond¬ 
ing to oscillatory Raman processes under the external 
drive, moreover, the resonant frequencies will be shifted 



Fig. 3: (Color online) Populations of photon numbers in 
dependence of the detuning. The parameter are: x /7 = 20; 
12/7 = 5, curves ( 1 ) and 17/7 = 20 , curves ( 2 ). 


due to the Stark effects as well. Indeed, in the second- 
order approximation of the perturbation theory for the 
interaction of an anharmonic oscillator with a monochro¬ 
matic field, the shift of oscillatory energy En in transi¬ 
tions through the states \n — l) and |n + 1 ) is obtained as 

AE — ( _ - _\ 

Va;+x(2n-l) a;+x(2n+l)F 

As the calculations based on Eqs. Q, and (§ show, 
the results strongly depend on all the parameters: the 
nonlinearity, the detuning and the amplitude of the driv¬ 
ing field. Below, we focus on comparative analysis of the 
regimes of strong and very strong nonlinearities in order 
to demonstrate selective excitations in KNR due to for¬ 
mation of well-resolved multi-photon resonance spectral 
lines. In this way, we concentrate on the particular cases 
with the values of nonlinearity: x /7 = 2 and x/x = 20 , 
considering the interplay of the ratio x/x^ Ihe detuning 
A/x and the amplitude of driving field iX/y. 


3.1 Non-selective excitations of oscillatory states 


In this subsection, we consider excitations of KNR on the 
base of Eqs. 0 and assuming the quantum regime, 
but not the formation of well-resolved energy levels of 
KNR. At first, two-dimensional plots of quantum mechan¬ 
ical mean photon numbers are depicted on Eig. in de¬ 
pendence of the detuning values A/x and the driving am¬ 
plitude iX/y. Comparing Eig. p!a|and Eig. lb we conclude 
that the maximum values of the mean photon numbers in 
the resonator are essentially small for large value of the 
ratio x/x = 20 in comparison with the case of x/x = 2 . 
This means that a high nonlinearity with respect to dissi¬ 
pation leads to formation of a strong quantum operational 
regime. 

Next, we analyze the photon number distribution in 
dependence of the detuning, for the case of a moderate 
nonlinearity, x/x = 2 . In this case the driving field is 
simultaneously resonant with many transitions between 
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Fig. 4: (Color online)Photon probability distributions Pq 
and Pi (a), P 2 and P 3 (b) in dependence of the nonlinear¬ 
ity. The parameters used are: Aj^ = —20, C/y = 20. 


consecutive oscillatory levels. The results of calculations 
the base of Eq. Q are depicted on Fig. 2a 


on 


as two- 

dimensional plots of the Fock state populations. These 
plots illustrate the formation of Poissonian-type distribu¬ 
tions from the thermal distribution, if the frequency of the 
driving field is tuned around the frequency of KNR. Nev¬ 
ertheless, these distributions are narrow compared to the 
Poissonian distributions for a higher amplitude of the driv¬ 
ing field. Photon-number distribution functions are plot¬ 
ted on Fig. 1^ for different values of the detuning. 


3.2 Selective resonant excitations of oscillatory states 

In this subsection we discuss the case of very strong non- 
linearities for which energy levels of KNR are well re¬ 
solved, and hence selective excitations of the oscillatory 
mode take place. Figure illustrates the dependence of 
Po, A, P 2 , P 3 populations of the photon-number states \n) 
on the detuning A/^ for the case of very strong nonlin¬ 
earity, xh — 20 . As we see, selective resonance excita¬ 
tions of the Fock states |1), |2), |3) take place for both 
cases of driving field: i ?/7 = 5 and Qjx = 20. Indeed, 
the populations Pi depicted on Fig. display a maxi¬ 
mum of 0.5 at A = 0, that corresponds to the resonant 
one-photon transition |0) ^ |1) at the frequency uj = ujq. 
The other maximum Pi = 0.38, for Q/x = 5, corresponds 
to the population of state |1) through the Raman process 
with energy conversation Eq + 2uj2 = cj/c + Pi . This pro¬ 
cess involves a two-photon excitation of the Fock state |2), 
(transition |0) ^ |2) at the frequency of the pump field 
co ’2 = cjo + x): and the decay |2) ^ |1) at the frequency 
= ^0 p 2 %. 


Fig. 5: (Color online)Photon probability distributions Pq 
and Pi (a), P 2 and P 3 (b) in dependence of the external 
force amplitude. The parameters are as follows: A/j = 

-20, x/7 = 20. 



Fig. 6 : Photon probability distributions Pq (a). Pi (b), P 2 
(c) and P 3 (d) in dependence on the detuning and the field 
amplitude. The nonlinearity parameter is: x /7 = 20. 


As we see on Fig.[^ for a comparatively weak driving 
field, 12/7 = 5, the population P 2 only reaches the value 
0.2 due to a two-photon resonant excitation | 0 ) ^ | 2 ) at 
the frequency Co ’2 = ^0 + X- the range of strong excita¬ 
tions the role of multi-photon processes increases. We ob¬ 
serve that these populations display peaked structures and 
the maximum values are realized for definite parameters of 
the detuning in accordance with resonant frequencies. In¬ 
deed, for the case of a strong field, C/y = 20, the Raman 
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processes of state |2) become effective. Particularly, the 
peaks at Aj^ = —40 and Aj^ = —60 correspond to the 
Raman processes through the states |3) at the frequency 
= ^0 + 2 % and |4) at co ’4 = cjo + 3% , respectively. As 
for weak excitation, 17/7 = ^ (I)? 


3d 


we see on Fig. 

population P 3 is negligible for all values of the detuning, 
while for the case of a strong excitation, 17/7 = 20 ( 2 ), a 
selective population takes place. 

In general, from analytically obtained results and the 
numerical analysis we can conclude that low-energetic states 
|1), 12) and |3) can be selectively excited with high effi¬ 
ciency in case of very strong nonlinearities. Nevertheless, 
we strongly demonstrate that values of Fock state popu¬ 
lations are limited in the cw steady-state regime of KNR. 
Indeed, on Fig. the population Pi is limited by the 
value 0.5, and the populations P 2 and P 3 are less than Pi. 
This effect of selective excitations can be interpreted as 
blockade of two or more photons in KNR by production of 
a single-photon Fock state. Thus, observation of these lim¬ 
its on the populations restricts the possibility to observe 
photon blockade effects in an over transient steady-state 
regime. 

We demonstrate this fact on Fig.[^ showing the depen¬ 
dence of populations from the parameters of nonlinearity 
at a fixed detuning value Aj^ = —20 and an external field 
amplitude 17/7 = 20. This regime involves two-photon res¬ 
onance excitation of KNR for x /7 = 20. It can be observed 
that effective excitations of the oscillatory mode take place 
only for particular values of the nonlinearity parameter, 
as a result the values of populations are limited, partic¬ 
ularly Pi is limited by 0.5, while the populations P 2 and 
P 3 are limited by the values 0.35 and 0.3 correspondingly. 
Analogously, the populations versus the amplitude of driv¬ 
ing field depicted on Fig. [^demonstrate limited behavior 
in steady-state regimes as well. However the behavior of 
populations of P 2 and P 3 depicted on Fig. |5b| differ from 
Po and Pi (Fig. 5a). The complete image of the popu¬ 


lations Po, Pi, P 2 and P 3 for the nonlinearity parameter 
x /7 = 20 (case of selective excitations), in dependence of 
the detuning A/j and the amplitude of the driving field 
17/7 are presented on Fig. 


3.3 Quantum statistics and Wigner functions 

In the end of this section, we turn to investigation of pho¬ 
ton statistics on the base of the normalized second-order 
correlation function Q for zero delay time. The corre¬ 
lation function describes the ratio of the simultaneously 
emitted photon pair number to the product of the in¬ 
dependently emitted photon number, and, on the other 
hand, is expressed through the variance of the photon 
number fluctuations as ((Z\n)^) = (n) + — 1 ). 

Thus, the condition < 1 corresponds to the sub- 
Poissonian statistics {{An)^) < {n) , the condition ^^(0) = 
1 corresponds to the Poissonian statistics and the condi¬ 
tion ^^(0) > 1 to the super-Poissonian statistics. The re¬ 
sults for g^ ( 0 ) and the mean photon numbers are depicted 
on Fig. for three values of the nonlinearity x /7 = 2, 




Fig. 7: The correlation function (a) and the mean exci¬ 
tation number (b) in dependence on the detuning for the 
amplitude of the driving field f2/j = 20. The nonlinearity 
parameter is: x/l = 2 (1), x/l = 10 (2), x/l = 20 (3). 


x /7 = 10 , x /7 = 20. As we see, for the regime of non¬ 
resonant excitations, (x/7 = 2), the correlation function 
approximately equals to unity for a large range of the 
detuning. This result is in accordance with the results 
depicted on Fig. and reflects the Poissonian statistics 
of KNR. Note, the correlation function displays a peak 
around the value of detuning Ajx = —35 in the critical 
or threshold region of photons production. An analogous 
result on critical growth of the correlation functions in the 
region of the generation threshold has been demonstrated 
for the non-degenerate parametric oscillator 43 . For the 
regimes of very strong nonlinearities the correlation func¬ 
tion is less than unity < 1), i.e. the mode displays 
sub-Poissonian non-classical statistics for a large range of 
detuning values from A/x = —40 to A/x = 40 due to 
selective excitations of KNR. 


At the end of this subsection we turn to discussion 
of phase-space properties of KNR in regimes of selective 
excitations. Typical results for the Wigner function of the 
cavity mode and the corresponding distributions of photon 
number states are depicted on Fig. and Fig. At first, 
we consider the case of the parameters f2/x = 5, A/x = 0, 
and x /7 = 20. This regime interests us in terms that only 
vacuum and single photon states are allowed here and the 
probabilities of both states equal to ( 0.5) (see. Fig. |^. 
The results of calculations (4) and (11) are illustrated on 
Figs. [8a| and [81^ As the Wigner function is positive over¬ 
all the phase-space, so KNR driven by a monochromatic 
force can not display any quantum-interference pattern. 
Nevertheless, the Wigner function displays non-classical 
photon-number effects in strong quantum regimes. On Fig. 
[^the Wigner function and the photon number probabili- 
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Fig. 8: The Wigner function (a) and the distribution of 
photon number states (b) for the parameters = 5, 

= 0 , xh = 20. 


Fig. 9: The Wigner function (a) and the distribution of 
photon number states (b) for the parameters 17/7 = 20, 
Z\/7 = -40, xh = 20. 


ties are illustrated for the case of a strong field 17/7 = 20 
and for the parameters Z \/7 = —40 and x /7 = 20. The 
probability of the vacuum state is neglectably low in this 
case (see, Fig. 3a). In this regime only one-photon and 


two-photon Fock states are effectively observed and the 
Wigner function displays a specific form in the phase- 
space. Thus, Fig.9 depicts distributions of n-photon states 
showing considerable contributions from the one- two- and 
three-photon states. The similar situation was already dis¬ 
cussed, although for the perturbation regime [^. 


4 Conclusion 

In conclusion, we have investigated various regimes of mode 
excitations in KNR driven by a continuous wave field for 
strong Kerr nonlinearities. The results have been obtained 
on the base of exact quantum solution of the Fokker- 
Planck equation in the complex P representation. This ap¬ 
proach allows non-perturbative consideration of KNR for 
various quantum operational ranges including cascaded 
multiphoton processes between oscillatory states stimu¬ 
lated by a strong driving field and for quantum noise of 
arbitrary strength. The selective excitations of the res¬ 
onator mode are realized in the cases of well resolved en¬ 
ergy spectra of an anharmonic oscillator by tuning the 


frequency of the driving field for wide ranges of the sys¬ 
tem parameters, the nonlinearity constant, the detuning 
and the amplitude of the driving field. On the base of this 
approach the mean photon numbers have been obtained 
as well as the photon statistics of the resonator mode has 
been analyzed using the probability distributions of pho¬ 
tons, the second-order photon correlation functions and 
the Wigner functions. It has been shown that for mod¬ 
erate nonlinearities photon state distributions are close 
to Poissonian distributions, while for strong nonlinearities 
selective excitations of photon number states may arise. 
Particularly, we have demonstrated the limited behavior of 
photon number-state populations in steady-state regimes 
due to concurrency of cascaded processes and decoherence 
effects. We observed that the population of one-photon 
state is limited by the value 0.5 while the population of 
two-photon state is limited by 0.35 for all ranges of the 
parameters. This conclusion restricts the possibilities for 
realization of the photon blockade due to cw excitations 
in over transient steady-state regimes. 
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8 


A R Shahinyan et aL: Selective excitations of a Kerr-nonlinear resonator: exactly solvable approach 


All authors contributed equally in the work presented 
in this paper. 


References 

1. A. Kubanek, A. Ourjoumtsev, 1. Schuster, M. Koch, P. W. 
H. Pinkse, K. Murr, and G. Rempe, Phys. Rev. Lett. 101, 
(2008) 203602. 

2. L. S. Bishop, J. M. Chow, J. Koch, A. A. Houck, M. H. 
Devoret, E. Thuneberg, S.M. Girvin, and R. J. Schoelkopf, 
Nature Phys. 5, (2009) 105. 

3. A. Imamoglu, H. Schmidt, G. Woods, and M. Deutsch, 
Phys. Rev. Lett. 79, (1997) 1467. 

4. K. M. Birnbaum, A. Boca, R. Miller, A.D. Boozer, T. E. 
Northup, and H. J. Kimble, Nature 436, (2005) 87. 

5. L. Tian and H. J. Carmichael, Phys. Rev. A 46, (1992) 
R6801. 

6. A. J. Hoffman, S. J. Srinivasan, S. Schmidt, L. Spietz, J. 
Aumentado, H. E. T’ureci, and A. A. Houck, Phys. Rev. 
Lett. 107, (2011) 053602. 

7. C. Lang, D. Bozyigit, C. Eichler, L. Steffen, J. M. Pink, A. 
A. Abdumalikov, M. Baur, S. Eilipp, M. P. da Silva, A. Blais, 
and A. Wallraff, Phys. Rev. Lett. 106, (2011) 243601. 

8. A. Earaon, 1. Eushman, D. England, N. Stoltz, P. Petroff, 
J. Vuckovic, Nat. Phys. 4, (2008) 859. 

9. P. Rabl, Phys. Rev. Lett. 107, (2011) 063601. 

10. A. Nunnenkamp, K. Bprkje, and S. M. Girvin, Phys. Rev. 
Lett. 107, (2011) 063602. 

11. Y. X. Liu, A. Miranowicz, Y. B. Gao, J. Bajer, C. P. Sun, 
and E. Nori, Phys. Rev. A 82, (2010) 032101. 

12. A. Verger, C. Ciuti, and 1. Carusotto, Phys. Rev. B 73, 
(2006) 193306. 

13. A. J. Hoffman, S. J. Srinivasan, S. Schmidt, L. Spietz, J. 
Aumentado, H. E. Tureci, and A. A. Houck, Phys. Rev. Lett. 
107, (2011) 053602. 

14. A. Earaon, A. Majumdar, and J. Vuckovic, Phys. Rev. A 
81, (2010) 033838. 

15. G. H. Hovsepyan, A. R. Shahinyan, and G. Yu. 
Kryuchkyan, Phys. Rev. A 90, (2014) 013839. 

16. G. H. Hovsepyan and G. Yu. Kryuchkyan, EPJD 69, 
(2015) 64. 

17. A. Miranowicz, M. Paprzycka, Y. X. Liu, J. Bajer and E. 
Nori, Phys. Rev. A 87, (2013) 023809. 

18. S. Eerretti and D. Gerace, Phys. Rev. B 85, (2012) 033303. 

19. T. C. H. Liew and V. Savona, Phys. Rev. Lett. 104, (2010) 
183601. 

20. H. J. Carmichael, Phys. Rev. Lett. 55, (1985) 2790. 

21. M. Bamba, A. Imamoglu, 1. Carusotto, and C. Ciuti, Phys. 
Rev. A 83, (2011) 021802(R). 

22. X.W. Xu and Y. J. Li, J. Phys. B 46, (2013) 035502. 

23. A. Majumdar, M. Bajcsy, A. Rundquist, and J. Vuckovic, 
Phys. Rev. Lett. 108, (2012) 183601. 

24. W. Zhang, Z. Y. Yu, Y. M. Liu, and Y. W. Peng, Phys. 
Rev. A 89, (2014) 043832. 

25. X. W. Xu and Y. Li, Phys. Rev. A 90, (2014) 033809. 

26. O. Kyriienko, 1. A. Shelykh, and T. C. H. Liew, Phys. Rev. 
A 90, (2014) 033807. 

27. S. Eerretti, V. Savona, and D. Gerace, New J. Phys. 15, 
(2013) 025012. 

28. H. Z. Shen, Y. H. Zhou, and X. X. Yi, Phys. Rev. A 91, 
(2015) 063808. 


29. H. Elayac and V. Savona, Phys. Rev. A 88, (2013) 033836. 

30. D. Gerace and V. Savona, Phys. Rev. A 89, (2014) 
031803(R). 

31. W. Leon-acuteski, Phys. Rev. A 55, (1997) 3874. 

32. G. Yu. Kryuchkyan, L. A. Manukyan and N. T. Muradyan, 
Optics Communications 190, (2001) 245. 

33. W. Leohski and R. Tanas, Phys. Rev. A 49, (1994) R20. 

34. W. Leohski and A. Kowalewska-Kudlaszyk, Prog. Opt. 56, 
(2011) 131. 

35. T. V. Gevorgyan, A. R. Shahinyan and G. Yu. 

Kryuchkyan, Phys. Rev. A 79, (2009) 053828. 

36. T. V. Gevorgyan, A. R. Shahinyan and G. Yu. 

Kryuchkyan, Phys. Rev. A 85, (2012) 053802. 

37. T. V. Gevorgyan, A. R. Shahinyan, Lock Yue Chew and 
G. Yu. Kryuchkyan, Phys. Rev. E 88, (2013) 022910. 

38. T. V. Gevorgyan and G. Yu. Kryuchkyan, J. Modern Op¬ 
tics 60, (2013) 860. 

39. P. D. Drummond and D. E. Walls, J. Phys. A 13, (1980) 
725. 

40. G. Yu. Kryuchkyan and K. V. Kheruntsyan, Opt. Comm. 
120, (1996) 132. 

41. K. V. Kheruntsyan, D. S. Krahmer, G. Yu. Kryuchkyan 
and K. G. Petrossian, Opt. Comm. 139, (1997) 157. 

42. K. V. Kherunstyan, J. Opt. B 1 , (1999) 225. 

43. G. Yu. Kryuchkyan, K. G. Petrosyan and K. V. 
Kheruntsyan, JETP Lett. 63 No 7, (1996) 502-506. 

44. J. Bourassa, E. Beaudoin, J. M. Gambetta and A. Blais, 
Phys. Rev. A 86, (2012) 013814. 

45. M. Boissonneault, A. C. Doherty, E. R. Ong, P. Bertet, D. 
Vion, D. Esteve, A. Blais, Phys. Rev. A 85, (2012) 022305. 

46. E. R. Ong, M. Boissonneault, E. Mallet, A. C. Doherty, 
A. Blais, D. Vion, D. Esteve and P. Bertet, Phys. Rev. Lett. 
110 , (2013) 047001. 



